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^ ! Abstract 
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. Within the framework of a dynamical and physically transparent model devel- 

!>• \ oped earlier, we study the time evolution of various quark flavours in the baryon-free 

region in ultrarelativistic heavy ion collisions. We show that even under optimistic 

■ conditions, the quark- gluon system fails to achieve chemical equilibrium. 
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i-C \ Quantum Chromodynamics (QCD), the standard theory of strong interactions, 

predicts that at very high temperature and/or density, the colourless hadronic mat- 
On. ter dissolves into its coloured constituents, the quarks and gluons, such that the 
(— I \ bulk properties of strongly interacting matter are governed by these degrees of free- 

■ dom. Such a locally colour-deconfined phase of matter is known as Quark Gluon 
Plasma (QGP). It is expected that the temperature and density achievable in ul- 
trarelativistic collisions of heavy ions is favourable for QGP formation, although 

' transiently. 

Many aspects of this transition, e.g'., the order of the phase transition, the value 
of the critical temperature, signals of the transition, thermodynamic equilibration 
, are still uncertain. So far, most of the works assume thermodynamic equilibrium 
throughout the evolution history, after an initial proper time Tj (~ 1 fm/c). Thus 
the thermodynamic quantities like pressure, entropy and temperature have direct 
meaning with respect to an equilibrated system. It has however been recently shown 
that the approach to equilibrium in ultrarelativistic heavy ion collisions proceeds 
through a succession of many time scales |jl], |2[. In particular, the dominance of the 
gg cross section over the qg or qq cross sections was argued to imply that the 
gluons equilibrate among themsleves appreciably earlier than the whole system of 
quarks and gluons. In a recent work [Q], referred to as 1, we showed that the time 
scale of kinetic equilibration for the light quarks (n, d) is ~ 1 fm/c; for massive 
quarks, it increases with the mass to rather large values. The significance of these 
considerations for QGP diagnostics has been discussed in 0, the emission of 
particles from the pre-equilibrium era (r^ < r < Tth), Tth denoting the time for full 
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thermodynamic equihbrium, may indeed populate the invariant mass or pT windows 
thought to be relevant for signals from the thermalised QGP. 

In this letter, we study the chemical evolution of the quarks in the system 
from the epoch Tg ( the proper time when gluons thermalise) onward. Chemical 
equilibration has very important implications for the signal of QGP formation |Q, 
^, ^. It would be ideal to study this problem in the QCD based kinetic theory and 
attempts along these lines have been made recently Q. Such calculations, appealing 
as they are, rely heavily on rather involved numerical simulations and also suffer 
from dependence on parameters needed to mimic the non-perturbative effects. We 
propose, as in 1, a physically transparent scenario which retains the essence of the 
kinetic theory to the following extent : the gluons thermalise completely at a proper 
time Tg earlier than the quarks. The gluons carry about half of the total momentum 
and the quarks carry a tiny fraction. We restrict ourselves to the central rapidity 
region where the number of valence quarks is assumed to be negligibly small. 

In order to treat the equilibration of the species one must follow the micro- 
scopic evolution of the phase space distribution function f{x^,p^), governed by the 
Boltzmann transport equation : 

L{f,} = C{f,} (1) 

where L is the Liouville operator and C the collision operator, the subscript q 
denoting a quark flavour. 

The Boltzmann equation is a non-linear, integro-differential equation for the 
phase space distribution function of the particles. For our present purpose, how- 
ever, the problem can be addressed in a rather simple manner without losing much 
information on the time scales for chemical equilibrium. 

In ref. ||8|, |9| the decoupling of the relic particles in the early universe was studied 
by integrating the Boltzmann equation over the momentum of the particle to obtain 
an equation for the evolution of number density. We use a similar approach to look 
at the number density evolution of the partons in ultrarelativistic collisions of heavy 
ions; the major difference between the two cases is that the expansion dynamics in 
the early universe {big bang) is governed by the Hubble expansion Q whereas in 
the case of heavy ion collision (mini bang), the expansion dynamics is assumed to 
be governed by the boost invariant (Bjorken) solution of realtivistic hydrodynamics 



[10|. The equation for the number density evolution then reads 



where Ug is non-equilibrium quark density, r the proper time and gg the statistical 
degeneracy of the quarks of flavour q. 

The most important contributions to the collision term in eq. (2) are from 



gluon fusion {gg <-> qq) and gluon decay {g — qq) [|lT|, |2|. Since we have assumed 
complete equilibration for the gluons, reactions like gg — > ggg... are not included in 
our work, the gluon density being determined by the temperature of the system at 
each instant. We ignore the quantum effects but take into account relativistic effects 
because of the very high temperature. Thus gluons are described by relativistic 
Maxwell - Boltzmann statistics, fg{E) = g/{2TT)^e~^^'^ . Then eq.(2) can be written 
as. 
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where neq{T) = gq/{2'K)^ j S-pjie^l'^ + 1), is the equihbrium density. Eq.(3) is a 
particular form of the Riccati equation; there is no closed form solution of this 
equation. 

In deriving eq.(3) we have assumed T (or CP) invariance, i.e., \M\ 



\M\ 



\qq—,gg- Rgg-^qq{T) IS the quark production rate per unit four volume by the 
reaction gg qq and Rg^qq gives the same for the gluon decay, g ^ qq [0, 12 1. The 
thermal quark masses have been taken into account [Q] using 77z?(T) = mq+g^T"^ /4.5 
in all the rates which was neglected in the work of Biro et al . The full derivation 
of the above equation has been omitted here for the sake of brevity; it will be 



reported elsewhere |14]. Eq.(3) is a very convenient form of the Boltzmann transport 
equation with all the relevant features for our present purpose. The first term on 
the right hand side of the equation represents the "dilution" of the density due 
to one dimensional expansion, the destruction of the quarks is propotional to the 
annihilation rates of qq and destruction is balanced by creation through the inverse 
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. The last term stands for the production of quarks due to 



gluon decay for which there is clearly no inverse process. 

To solve eq.(3) for ng(r) we need to know r(r), which is obtained by solving 
the energy conservation equation 



de e + P 
dr T 
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where P is the pressure and e the energy density. An equation of state relating P 
and e can be parametrized in the form 



P 



Then, for the non-equilibrium configuration ||T5|, we can write 



e = [arg + 6r,(T)] 



(5) 
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where a = (87r^)/45c7^ and b = {7tt'^) / 60 Nf Cg"^ . In chemical equilibrium, rq{= 
iT'q/^eq) for all spccics should be unity and in correspondence with our premise, we 
take rg =1 all through. A'^^ is the effective number of massless flavours; appearance 



of quark masses amounts to reducing the actual number of flavours |16, 17]. For 
the present, we confine our attention only to the three lightest flavours, u,d and s 

We then have. 
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To solve eq.(3), we need to specify the initial conditions. The initial values of 



nq{Tg)^s are obtained from the integration of the structure function [19| 
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where A is the mass number of the cohiding nuclei (208 in our case for Pb), 
is the minimum value of x, taken to be 0.02 ^, dictated by the applicability of 
perturbative QCD. The values of ng{Tg) = 1.42,2.54 and 3.05 for SPS, RHIC and 
LHC respectively. The corresponding initial values of the temperature (Tg) are 330, 
500 and 660 MeV at proper times rg=0.54, 0.3 and 0.25 fm/c, respectively 0. 

The solutions of eqs.(3) and (7), obtained by the Runge-Kutta method, provide 
us with the non-equilibrium density rig^r) and the cooling law, 



a + brq{Tg) 
a + brq{T) 



n3/4 

T^rg (9) 



As is evident from the above discussions, the non-equilibrium density Ug has an 
explicit dependence on r and an implicit dependence on r through T{t). But the 
equilibrium density neg{T) has only an implicit dependence on r through T[t). 
The ratio thus assumes an universal feature, since the implicit time dependence 
gets eliminated. The time dependence of the ratio Vq can then be used as a ready 
marker for chemical equilibrium; the time at which the explicit time dependence of 
Tq vanishes corresponds to the time for chemical equilibration for the flavour q. 

In contrast to the earlier work of Biro et al [^] , we have not included reactions 
like qq QQ, since the initial system is dominated by gluons[Q, ^; the quark density 
is very low compared to the gluons. These reactions are suppressed by a factor of 
about 1/16 M compared to the case when the quarks are in complete equilibrium. 



Biro et al |1S] also estimated their initial density from the HIJING model 

The ratios = Uq/ueq are plotted in fig.l for SPS, RHIC and LHC energies. 
At RHIC and LHC energies, the ratios for u{d) quarks and s quarks are the same at 
early times but at late times the u{d) quarks dominate over the s quarks. The reason 
is that at early times the thermal mass (oc T) dominates over the current quark mass 
but at later times the thermal mass becomes small and the main contribution comes 
from the current quark mass which is large for s quarks (m^ ~ 150Mey) compared 
to the u{d) quark masses {rriy^i^d) ~ l^MeV); there occurs thus some Boltzmann 
suppression for s quarks relative to u{d) quarks. At SPS energies, at early times 
the ratio for s quarks is large compared to the u{d) quarks, this is due to the initial 
normalisation. The non-equilibrium density n^(^) ~ Ug but the equilibrium density 
for s quarks is small compared to u and d quarks due to larger mass. 

At small time scales (larger temperature) the slope of the ratio for LHC energies 
is opposite to that of SPS and RHIC. The non-equilibrium density contains the effect 
of expansion as well as production, where as the equilibrium density is affected by 
expansion only. In the ratio, the effect of production dominates over the expansion 
at LHC i.e. larger temperature, so the ratio increases. However at RHIC and 
especially at SPS the quark production is not so intense to dominate over expansion. 

In fig.l, we have also shown the time variation of temperature, given by the 
cooling law (eq.(9). The lifetime of the quark-gluon system is the time when T falls 
below Tc, taken to be 160 MeV as an optimistic estimate. It is abundantly clear that 
under either of the two criteria, ^ 1 and/or flatness of - the quarks remain out 
of chemical equilibrium in the baryon-free region at SPS, RHIC and LHC energies. 
Detailed account of these effects in QGP diagnostics must be considered; work along 
this line is in progress. 
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Figure Captions 

Figure 1: Ratio of non-equilibrium density to equilibrium density and temperature 
as functions of time at SPS, RHIC and LHC energies. 
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